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Abstract 

We study the effect of dissipation on the density profile of a one-dimensonal 
gas that is subject to gravity. The gas is in thermal equilibrium at temperature 
T with a heat reservoir at the bottom wall. Perturbative analysis of the 
Boltzmann equation reveals that the correction due to dissipation resulting 
from inelastic collisions is positive for < z < Zc and negative for z > Zc 
with z the vertical coordinate. The numerically determined value for Zc is 
rngZc/ksT ~ 1.1613, where g is the gravitational constant and m is the 
particle mass. 

I. INTRODUCTION 

Granular materials are basically a collection of meso- to macroscopic particles that in- 
teract with each other via short range repulsive potentials [1]. For this reason, they may be 
considered as a very dense molecular gas, but with several key differences: first, since the 
particles are macroscopic, gravity plays an important role and cannot be dismissed, as it 



*E-mail address: kurtze@golem.phys.ndsu.NoDak.edu 
^E-mail address: dh09@lehigh.edu 

1 



can for a molecular gas. In fact, most of the unique features of granular materials disappear 
in the absence of gravity. Second, collisions among the particles are inelastic, so dissipation 
must be included in studies of the response of granular materials to external stimuli. It 
has recently been observed that dissipation appears to be responsible for the existence of a 
uniform cooling state and the associated clustering instability [2-7], and modifies the veloc- 
ity distribution function from the usual Maxwellian to a power law or exponential function 
[8,9]. The question of whether the appearance of the clustering instabihty is a generic phe- 
nomenon or an artifact of the hard sphere potential remains open; nonetheless we note that 
the role of dissipation in nonequilibrium dynamics has always been subtle, complex, and 
mathematically challenging. 

In this paper, we examine the effect of dissipation on the density profile of a one- 
dimensional granular gas within the framework of the Boltzmann equation. In so doing, 
we implicitly treat the grains as point particles, and thus ignore one of the crucial aspects 
of granular materials, namely the excluded volume effect [10,11]. Neverthless, the outcome 
of this investigation is interesting and requires further study because it makes a nontrivial 
prediction: the dissipation leads to an increase in density near the source and a decrease 
away from the source, with the dividing line appearing at a dimensionless height gz/T of 
approximately 1.1613. This may have some relevance to the recent observation of a cluster- 
ing instability near the source in the presence of gravity [7] . However, it remains to be seen 
whether our result persists in higher dimensions, or survives the inclusion of the excluded 
volume effect. 

II. PERTURBATIVE ANALYSIS OF THE ONE-DIMENSIONAL BOLTZMANN 

EQUATION 

Consider a one-dimensional gas of N identical classical point particles moving along the 
z axis above a wall at z = 0, and acted upon by a constant gravitational field g in the 
negative z direction. The particles collide with one another inelastically, with a coefficient 
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of restitution r < 1. Let Nf{u, z) dudz be the average number of particles found between 
heights z and z+dz with velocity between u and u+du, in the steady state. This distribution 
function is given by the one-dimensional Boltzmann equation with dissipation [2], 

f d d\ 

\dz ~ ^du J " ^' ~ 

The collision terms on the right hand side are the "creation" term, 

C{u, z;q)^ J J \u' - u"\ f{u', z) f{u", z) 5[u - qu' - (1 - q)u"] du' du", (2) 

which gives the rate at which particles collide at z and come out of the collision with velocity 
u, and the "annihilation" term, 

A{u, z) = f{u, z) j \u' — u\ f{u' , z) du' = C{u, z;0), (3) 

which is the rate at which particles with velocity u collide at z and emerge with some other 
velocity. The parameter q is related to the coefficient of restitution by r = 1 — 2g. 

We assume that the bottom wall at 2; = maintains the particles there in some given 
velocity distribution fb{u). This roughly means that the particle distribution / obeys the 
boundary condition f{u, z = 0) = fb{u) at the bottom wall, but we must be a bit more 
careful than this. Since the particles are not allowed into the region of negative z, the 
bottom wall can only propel particles upward. Particles ai z — 0+ with positive velocities u 
must have been launched by the bottom wall, but those with negative velocities must have 
either undergone collisions or at least been in flight for a finite time since they were last in 
contact with the bottom wall. Therefore the velocity distribution of these latter particles is 
not fb{u), but instead is something which was established by the dynamics described by the 
Boltzmann equation itself. Only the former particles, then, those with positive velocities, 
must be in the velocity distribution imposed by the bottom wall. Thus the correct boundary 
condition is 

f{u,0)^Mu) for u>0. (4) 

The other boundary conditions are the obvious ones, namely that / should vanish for 2; 00 
and for \u\ — > 00. 



We now define a new parameter 
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e = 



(5) 



1-q 1 + r' 

which is small when the collisions are only slightly dissipative, i.e., when the coefficient of 
restitution is only slightly less than unity. Writing the creation term (2) in terms of e and 
carrying out the integral over u" gives 

C{u, z; g) = (1 + j \u' - u\ f{u', z) f{u + e[u - u'],z) du'. (6) 

We now expand this in powers of e and substitute it into the Boltzmann equation (1) to get 

du'+ 



2fiu,z) + iu-u')^ 



\u' — u\ f{u', z) 



du' ^ — . 



(7) 



This is the starting point of our calculation. 

Let us write the solution of (7) as a power series, 

/(m, z) = /o(m, z) + e/i(M, z) + 6^2 (m, z) + 



(8) 



At zeroth order, /o must satisfy 



dz du 



This is satisfied identically by any function of the energy variable 

E{u, z) = ^u'^ + gz. 



(9) 
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Thus to zeroth order we may satisfy the Boltmann equation and the boundary condition (5) 
by writing 

fo{u,z) = M^2E{u,z)). (11) 

This is the complete solution to the problem in the case of elastic collisions. Note that there 
would have been no solution to this problem had wc attempted to impose the boundary 
condition (5) for all velocities u, with a function which is not even. 



We will most often be interested in an assemblage which is in thermal equilibrium at 
z — 0. In this case we have 

fo(u,z) ^XeM-E/T) =A^exp(-(ii2 + 2yz)/2T), (12) 

where the temperature is mT/kB, where m is the particle mass and is Boltmann's con- 
stant, and the normalization constant jV is 

= g/V2^. (13) 

The zeroth order density profile is obtained by integrating over velocity, 

Po(z)^N fo{u,z)du. (14) 

J CO 

In the thermal equilibrium case (12), this becomes 

Po{z) = {Ng/T)exp{-gz/T). (15) 

We now calculate the first-order correction fi to the distribution function. Prom (7), we 
see that this must satisfy 

r r ;^f„l 

du'. (16) 



2f,iu,z) + iu-u')^ 



The right side of this equation is the li-derivative of fo{u, z) J {u — u')\u — u'\fo{u', z) du', but 
we have not found this intriguing fact to be very useful for our purposes. Note, however, 
that fo{u) is an even function of u, and therefore the right side of (16) is also even in u. 
Since the operator on the left side is odd in u, this implies that solving (16) will lead to 
a function which is odd in u. Since we find the density profile by integrating over all u, 
this function will not change the density from the zeroth order result (14). However, the 
operator in (16) also annihilates all functions of E{u, z), so /i can also contain any function 
of E. Such a function is even in u, and so does contribute to the density profile. 

For the thermal equilibrium case, we can calculate /i explicitly. First note that the z 
dependence of the right side of (16) consists entirely of a factor exp{—2gz/T) coming from 
the two /o factors in the integrand. To take advantage of this, we write fi in the form 

h{u, z) = hi{E{u, z)) - g-'f^{E{u, z))F,{u). (17) 



Note that Fi is independent of z. Substituting this expression, and the exphcit form (12) of 
/o, into (16) then gives 



du 



/CO 
W- 
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2 + 
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u{u' — u) I U" — U 



expl (18) 



T 

Changing integration variables from u' to v = u' — u and doing some algebra, we find that 
this can be written as 

AT = - (4 --jexp(^--j cosh (19) 

SO we finally have 

F,{u) = 2T f (4 - "pj exp (-^ j sinh (|^) dv. (20) 

This can also be written out explicitly in terms of error functions. 

We must now determine the function hi{E), which forms the even part of fi. This is 
slightly more subtle than it first appears. The first thing we should do is apply the boundary 
condition (4) at 2; = 0, thus maintaining the thermal equilibrium distribution of particles 
coming upward from the bottom wall. Thus we would choose hi{E) to be g~^fQ{E)Fi{\/2E). 
However, as we will see below, this has a problem: the resulting density profile has fewer 
than A'^ particles. In order to remedy this situation, we must change the normalization of 
/o, or equivalently include an extra term proportional to /o(-E) in hi{E), with a coefficient 
chosen to make the total number of particles in the system again equal to N. This amounts 
to replacing the boundary condition (4) with a slightly more general condition of the form 

f{u,0)(xMu) for u>0, (21) 

with the proportionality constant chosen so that / is properly normalized. Thus we maintain 
the correct number of particles, and the functional form of the distribution of particles 
launched from the bottom wall. This relaxation of condition (4) is physically reasonable, 
because the normalization of the distribution is not directly observable. 

We now have an expression for the distribution function / correct to first order in the 
dissipation parameter e, 

f{u, z)^{l + eN,)fo{E) + eg-^fl{E)[F,{^) - F,{u)] + 0{^), (22) 
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where E is given by (10), /o by (12) and (13), Fi by (20), and the constant jVi is yet to 
be determined. The density profile is obtained by integrating over u; since Fi{u) is an odd 
function it drops out, leaving us with 

N-'p{z) = (1 + eJ\fi){g/T) exp{-gz/T) + eg-' J f^{E)F,{V2E) du + 0{e^). (23) 

Integrating this equation over all z gives unity on the left side, so we see that Afi must be 
given by 

M = -g-' / dz / dufl{E)Fi{V2E). (24) 

Jo J-oo 

Remarkably, the integrals on the right can be evaluated analytically, giving the result 

Ml = 2/n. (25) 

Using this result in (23), we obtain the central result of this paper: an expression for the 
first-order correction to the density profile due to dissipation. After rescaling the integration 
variables u and v by factors of T^/^, we get 

TT I 

roo roo , "| 

+ exp(-2c/^/r) du dv{4 — v'^) exp{—u^) exp(— /2)smh{vJu'^ + 2{gz/T))\ . (26) 

We have evaluated (26) numerically as a function oi gz/T after first doing the integral 
over V analytically (in terms of error functions). The result is shown in Figure 1. The most 
important feature of pi{z) is the fact that it is positive near the source and negative away 
from the source, with a dividing line at Zf. = 1.1613T/5f. This is somewhat larger than T/g, 
which is the average z coordinate of the particles when their collisions are elastic. Thus 
below Zc the dissipation enhances the density profile, while above Zc it suppresses it. We 
might expect this behavior, since the dissipation in the collisions removes mechanical energy 
from the system. The effect of the thermal reservoir at the bottom wall is independent 
of the characteristics of the coUisions, so the particles leaving the bottom wall are in the 
same distribution whether the collisions are elastic or not. However, when their collisions 
are inelastic, every collision reduces the total mechanical energy of the assemblage, thus 



reducing the average height which the particles reach. This may have some relevance to 
the recent experimental observation of the clustering instability by KudroUi et al. [7] , who 
noticed the migration of the clustering instability along the direction of gravity and toward 
the driving source. 

We note also that the perturbative correction to the density profile reaches a minimum 
of —0.206 at gz/T ^ 2.22, and it also has a small maximum of 1.0261 at gz/T ^ 0.037. We 
have not yet identified the physical mechanism responsible for the presence of this maximum. 

It remains to be seen whether the enhancement of the density profile near the source in 
the presence of dissipation will persist in higher dimensions, or if it will survive the inclusion 
of excluded volume effects in the model. Future studies that elucidate these questions may 
provide interesting new insights into the role of dissipation in granular dynamics. 
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Figure Captions 

Fig. 1. The first-order correction pi{z) to the density profile, given by (26), plotted as 
a function of the dimensionless height z = gz/T. The value at 5 = is 1. The graph then 
rises to a small maximum z — 0.0372 with a height of 1.0261, crosses zero z — 1.1613, 
and reaches a minimum aX z — 2.2217 with a value of —0.2057. For large z it approaches 
the axis as — 25exp(— 5). 
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